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ABSTRACT Bacterial DNA and live bacteria have been detected in human urine in the absence of clinical infection, challenging the 
prevailing dogma that urine is normally sterile. Urgency urinary incontinence (UUI) is a poorly understood urinary condition 
characterized by symptoms that overlap urinary infection, including urinary urgency and increased frequency with urinary in- 
continence. The recent discovery of the urinary microbiome warrants investigation into whether bacteria contribute to UUI. In 
this study, we used 16S rRNA gene sequencing to classify bacterial DNA and expanded quantitative urine culture (EQUC) tech- 
niques to isolate live bacteria in urine collected by using a transurethral catheter from women with UUI and, in comparison, a 
cohort without UUI. For these cohorts, we demonstrated that the UUI and non-UUI urinary microbiomes differ by group based 
on both sequence and culture evidences. Compared to the non-UUI microbiome, sequencing experiments revealed that the UUI 
microbiome was composed of increased Gardnerella and decreased Lactobacillus. Nine genera (Actinobaculum, Actinomyces, 
Aerococcus, Arthrobacter, Corynebacterium, Gardnerella, Oligella, Staphylococcus, and Streptococcus) were more frequently cul- 
tured from the UUI cohort. Although Lactobacillus was isolated from both cohorts, distinctions existed at the species level, with 
Lactobacillus gasseri detected more frequently in the UUI cohort and Lactobacillus crispatus most frequently detected in con- 
trols. Combined, these data suggest that potentially important differences exist in the urinary microbiomes of women with and 
without UUI, which have strong implications in prevention, diagnosis, or treatment of UUI. 

IMPORTANCE New evidence indicates that the human urinary tract contains microbial communities; however, the role of these 
communities in urinary health remains to be elucidated. Urgency urinary incontinence (UUI) is a highly prevalent yet poorly 
understood urinary condition characterized by urgency, frequency, and urinary incontinence. Given the significant overlap of 
UUI symptoms with those of urinary tract infections, it is possible that UUI may have a microbial component. We compared the 
urinary microbiomes of women affected by UUI to those of a comparison group without UUI, using both high-throughput se- 
quencing and extended culture techniques. We identified statistically significant differences in the frequency and abundance of 
bacteria present. These differences suggest a potential role for the urinary microbiome in female urinary health. 



Received 15 May 201 4 Accepted 5 June 2014 Published 8 July 2014 

Citation Pearce MM, Hilt EE, Rosenfeld AB, Zilliox MJ, Thomas-White K, FokC, Kliethermes S, Schreckenberger PC, Brubaker L, Gai X, Wolfe AJ. 201 4. The female urinary 
microbiome: a comparison of women with and without urgency urinary incontinence. mBio 5(4):e0 1283-1 4. doi:10.1 1 28/mBio.01283-14. 
Editor Martin Blaser, New York University 

Copyright © 2014 Pearce eta I. This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported 
license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited. 
Address correspondence to Alan J. Wolfe, awolfe@luc.edu, or Xiaowu Gai, xiaowu_gai@meei.harvard.edu. 
This article is a direct contribution from a Fellow of the American Academy of Microbiology. 



The female urinary tract is a previously overlooked microbial 
niche. Recent detection of bacterial DNA (1-3) and live bac- 
teria (4, 5) in urine from women with and without lower urinary 
tract symptoms has provided compelling evidence that the female 
urinary tract possesses its own unique microbiota. Research into 
the roles of these bacterial communities in urinary health and 
disease requires sensitive and specific detection and classification. 
Although the comparatively low bacterial load in urine challenges 
such effort, expanded culture conditions have allowed researchers 
to detect bacteria at lower levels than previously used techniques. 
For example, we recently reported that bacteria overlooked by 



standard culture conditions are cultivable by adjusting the growth 
conditions to include increased urine volume, diverse growth me- 
dia and atmospheric conditions, and lengthened incubation time 
(4). Using expanded quantitative urine culture (EQUC), we iso- 
lated bacteria from 80% of examined urine samples collected by 
transurethral catheter, most (92%) of which had been reported as 
"no growth," using a standard clinical microbiology urine culture 
protocol and a 10 3 -CFU/ml threshold (4). Despite this significant 
advance, many urinary bacteria cannot be cultured under these 
conditions. Thus, culture-independent methods for bacterial de- 
tection, such as high-throughput sequencing of the 16S rRNA 
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TABLE 1 Participant demographics and symptoms 

Value for group («")'' or statistic 

Non-UUI Control 

Characteristic UUI (60) (58) Pvalue c 



Demographics 



A Pf ( VT<i) 


63 (±12) 


49 ("+14) 


<0.05 


Race [no. (%)] 






0.42 


Caucasian 


50 (83) 


41 (77) 




African- American 


10 (17) 


11 (21) 




Asian 


0(0) 


1(2) 




Estrogen status [no. (%)] 






<0.05 


Estrogen positive 


16 (27) 


31 (61) 




Estrogen negative 


43 (73) 


20 (39) 




Body mass index (kg/m 2 ) 


32 (±8) 


28 (±6) 


0.01 


Hypertension [no. (%)] 


21 (35) 


9(17) 


0.03 


Coronary artery disease [no. (%)] 


9(15) 


1(2) 


0.02 


Symptoms 








Symptom score'' 


66 (±23) 


4 (±7) 


<0.05 


Health-related quality of life (HRQL)'' 


56 (±25) 


99 (±3) 


<0.05 


Urinary Distress Inventory (UDI) e 


102 (±47) 


11 (±H) 


<0.05 


Pelvic Organ Prolapse Distress Inventory (POPDI) c 


70 (±52) 


28 (±32) 


<0.05 


Colorectal-Anal Distress Inventory (CRADI) C 


71 (±61) 


43 (±62) 


0.02 



" n, no. of subjects. 
6 Mean±SD or no. (%). 

c Pearson's chi-square and Fisher's exact tests were used with categorical variables. Student's f test was used with continuous variables. 

a Based on OAB questionnaire. 

e Based on Pelvic Floor Disease Inventory. 



gene, have emerged as the predominant research technique, espe- 
cially as they become increasingly accessible due to declining se- 
quencing cost and improved bioinformatics tools (6). 

There is a clinical need for these improved research techniques. 
For example, urgency urinary incontinence (UUI) is poorly un- 
derstood chronic urinary condition that is often attributed to ab- 
normal neuromuscular signaling and/or functioning. However, 
alternative or complementary mechanisms beyond neuromuscu- 
lar abnormalities must be considered, since UUI pharmacological 
treatments aimed at the overactive bladder are ineffective in ap- 
proximately half of the pharmacologically treated UUI population 
(7,8). 

UUI is common, affecting 1.5 to 22% of the U.S. population, 
most frequently women and the elderly (8). Affected patients ex- 
perience a sudden, intense need to urinate with involuntary urine 
loss that detracts from the patient's quality of life (9). A propor- 
tion of UUI cases resolve over time; however, it is a chronic con- 
dition for most women (8). There is a large socioeconomic bur- 
den, with the costs of UUI in the United States projected to reach 
$76.2 billion by 2015 (10). 

The clinical diagnosis of UUI requires exclusion of urinary 
tract infections (UTI); thus, infectious etiology is not considered 
for UUI. Given the clinical similarly of UUI and UTI symptoms, 
however, we and others have recently used expanded culture tech- 
niques (4, 5) and quantitative PCR (2) to show the presence of 
bacteria in standard culture-"negative" urine samples collected 
from UUI patients. 

In this analysis, we utilized both 16S rRNA gene sequencing 
and EQUC to characterize the microbiome in urine obtained by 
transurethral catheter from women seeking treatment for UUI 
and a comparison group of women without UUI. Utilizing both 
techniques, we identified statistically significant differences in the 



frequency and abundance of bacteria present, suggesting a poten- 
tial role for the urinary microbiome in female urinary health. 

RESULTS 

Cohort description. Table 1 displays the UUI symptoms by co- 
hort. Consistent with the intended composition of the cohorts, the 
UUI cohort reported more distress on both the urinary and pro- 
lapse subscales of the Pelvic Floor Distress Inventory (PFDI) (uri- 
nary distress inventory [UDI], 102 [±47] versus 11 [±11] [P < 
0.05]; Pelvic Organ Prolapse Distress Inventory [POPDI], 70 
[±52] versus 28 [±32] [P < 0.05] ). The UUI cohort also reported 
more-severe overactive bladder (OAB) symptoms (66 [±23] ver- 
sus 4 [±7];P< 0.05) and decreased quality of life as measured by 
health-related quality of life (HRQL) (56 [±25] versus 99 [±3]; 
P < 0.05) compared to the non-UUI cohort. 

Table 1 also displays the participant demographics by cohort. 
Marital status, diabetes, smoking, and prior pelvic surgeries did 
not differ by group. However, the UUI cohort was heavier (body 
mass index [BMI],32[±8] versus 28[±6];P = 0.01), less likely to 
be using estrogen (27% versus 61%; P < 0.05), and older (63 
[±12] versus 49 [±14]; P < 0.05). 

Sequence-based characterization of female urinary micro- 
biome. Bacterial DNA was detected in similar proportions of 
urine samples (UUI, 63.9% [23/36]; non-UUI, 65.8% [25/38]). 
The lack of bacterial detection in the remaining samples may be 
due to low bacterial load, insufficient bacterial lysis, and/or primer 
bias rather than being conclusive evidence of no bacteria in these 
urine samples. 

Further analyses were performed on the samples with detect- 
able bacterial DNA (for UUI, n = 23; for non-UUI, n = 25). The 
sequences from these samples were classified into 22 phyla, 34 
classes, 69 orders, 150 families, and 386 genera. The most abun- 
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FIG 1 Urinary microbiome profile by cohort based on 16S rRNA gene V4 sequencing. Stacked bar plots depict the sequence abundances of the 15 most 
abundant genus- or family-level taxa in the UUI and non-UUI cohorts. Taxa were ranked according to mean abundance across all samples. The y axis represents 
the percentage of sequences for a particular bacterial taxa; the x axis represents the study participants separated by cohort. The family Enterobacteriaceae could 
not be classified to the genus level. The remainder of sequences were combined in the category labeled Other. 
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dant phyla detected were Firmicutes, with a median abundance of 
60%, followed by Actinobacteria (16%), Proteobacteria (1.5%), 
and Bacteroidetes (0.6%). From each sample, the majority of se- 
quences could be classified to the genus level. The two exceptions 
(C031 and C033) were predominantly Enterobacteriaceae (Fig. 1). 
In the remainder of the samples, the percentage of sequences per 
sample that could not be classified to the genus level ranged from 
0.02 to 24.1%. At the genus level, the majority of urine samples 
were dominated by one or two bacterial families or genera, most 
frequently Lactobacillus and Gardnerella (Fig. 1). Urine samples 
with the same dominant taxa clustered together in a dendrogram 
generated via hierarchical clustering of the Euclidean distance be- 
tween urine samples (Fig. 2). We observed six groups, now termed 
urotypes, which were named based on either the dominant family 
or genus (Gardnerella, Sneathia, Staphylococcus, Enterobacteria- 
cae, and Lactobacillus) or the lack of a dominant family or genus 
("diverse"). The most frequent urotype in both cohorts was Lac- 
tobacillus (43% UUI and 60% non-UUI controls), followed by 
Gardnerella (26% UUI and 12% non-UUI controls) (Table 2). 
Three urotypes were present in both cohorts {Gardnerella, Lacto- 
bacillus, and diverse). The Staphylococcus and Sneathia urotypes 
were present only in the UUI cohort, whereas the Enterobacteri- 
aceae urotype was present only in the non-UUI cohort; however, 
these differences were not statistically significant at our current 
sample size (Fig. 1 and Table 2). 

We calculated the frequency that each genus was observed in 
each cohort (Fig. 3A). While Lactobacillus was detected in every 
sample in both cohorts, some genera were observed more fre- 
quently in one cohort than in the other. For example, Gardnerella 
(99% UUI, 60% non-UUI; P = 0.003) and Aerococcus (74% UUI, 
28% non-UUI; P = 0.002) were detected more frequently in the 



UUI cohort than in the non-UUI one. Conversely, Staphylococcus 
(61% UUI, 92% non-UUI; P = 0.01) was observed less frequently 
in the UUI cohort than among non-UUI women. Actinobaculum 
(61% UUI, 44% non-UUI; P = 0.24) and Sneathia (22% UUI, 4% 
non-UUI; P = 0.09) exhibited a trend toward increased detection 
in the UUI cohort; however, this difference was not statistically 
significant at our current sample size. 

By calculating frequency, all samples are weighted similarly 
regardless of sequence abundance. For example, urine samples 
collected from control participants C013 and C036 both con- 
tained Lactobacillus sequences; however, the sequence abun- 
dances differed. For C013, Lactobacillus represented more than 
99% of the sequences; for C036, Lactobacillus represented less than 
5% of the total sequences (Fig. 1). To further investigate these 
differences, the median sequence abundance was calculated and 
compared between the cohorts (Fig. 3B). The UUI cohort exhib- 
ited decreased Lactobacillus (14% UUI, 67% non-UUI; P = 0.01) 
and increased Gardnerella (4% UUI, 0% non-UUI; P = 0.003) 
sequence abundances compared to those for the control cohort. 
Figure S2 in the supplemental material displays the distribution of 
the 15 most abundant taxa detected by sequencing. 

To measure the richness of the UUI and non-UUI urinary 
microbiomes, the number of observed operational taxonomic 
units (OTUs) and the Chaol estimator were calculated. To mea- 
sure diversity, the Shannon index and inverse Simpson index were 
calculated. Overall, the urinary microbiome exhibited an average 
of 97 OTUs per urine, a Chaol index of 7,210, an inverse Simp- 
son's index of 1.72, and a Shannon index of 3.93 (Table 3). There 
were no statistically significant differences between the UUI and 
non-UUI control cohorts based on these estimators of richness 
and diversity. 
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FIG 2 Clustering of the urinary microbiome into urotypes. The dendrogram was based on hierarchical clustering of the Euclidean distance between samples in 
the combined UUI and non-UUI cohorts. The dashed line depicts where the clades were divided into 6 urotypes: Gardnerella, Sneathia, Diverse, Staphylococcus, 
Enterobacteriaceae, and Lactobacillus. The stacked bar plot below the dendrogram depicts the sequence abundances of the overall most abundant taxa. 



Culture-based characterization of the female urinary micro- 
biota. Of the urine specimens assayed via EQUC, 71/90 (78.9%) 
grew bacterial species, while 64 of the 71 (90.1%) that grew bac- 
terial species were deemed culture negative (no growth) by stan- 



dard culture. For the UUI and control cohorts, standard culture 
had false-negative rates of 90.3% and 90.0%, respectively. This 
highlights the limitations of the standard clinical microbiology 
protocol. 
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TABLE 2 Urotype frequency by cohort 

No. (%) of participants for group («") 







Non-UUI control 




Urotype 


UUI (23)'' 


(25) fc 


P value' 


Gardnerella 


6 (26) 


3(12) 


0.23 


Sneathia 


2(9) 


0(0) 


0.22 


Diverse 


4(17) 


5(20) 


0.99 


Staphylococcus 


1(4) 


0(0) 


0.48 


Enterobacteriaceae 


0(0) 


2(8) 


0.49 


Lactobacillus 


10 (43) 


15 (60) 


0.39 



" n, no. of subjects. 

b Analysis performed on the subset of samples with detectable bacterial DNA. 
c Fisher's exact test. 



Rarefaction curves comparing the number of bacteria species 
to the number of urine samples assayed reveal that the UUI cohort 
curve began to plateau at a higher number of isolates than the 
non-UUI cohort curve (Fig. 4), indicating that these two cohorts 
possess different urinary microbiota. In support of this supposi- 
tion, there was a statistically significant difference between the 
median numbers of bacterial isolates cultured from the UUI 
women and those for the non-UUI women (4 [interquantile ratio 
{IQR} = 1 to 7] versus 1 [IQR = 1 to 7]; P < 0.001). 

The frequencies at which we cultured each genus were com- 
pared between cohorts. Nine genera {Actitwbaculum, Actinomy- 
ces, Aerococcus, Arthrobacter, Corynebacterium, Gardnerella, Oli- 
gella, Staphylococcus, and Streptococcus) were more frequently 
isolated from the UUI cohort than from controls (Fig. 5). Four of 




Taxa 


UUI(N=23) 


Non-UUI Control (N=25) 


p-value 


Median (%) 


IQR (%) 


Median (%) 


IQR (%) 


Alloscardovia 


0 


0 


0 


0 


0.92 


Anaerococcus 


0.2 


(0-5.7) 


0.5 


(0-0.7) 


0.66 


Atopobium 


0 


(0-0.02) 


0 


(0-0.1) 


0.24 


Bifidobacterium 


0 


(0-0.1) 


0 


(0-0.03) 


0.65 


Corynebacterium 


0.1 


(0-0.7) 


0 


(0-3.8) 


0.09 


Enterobacteriaceae 


0 


(0-0.2) 


0.1 


(0-0.3) 


0.09 


Gardnerella 


4.0 


(0-50.1) 


0 


(0-0.1) 


0.003* 


Lactobacillus 


13.6 


(0-61.6) 


66.6 


(2.1-91.1) 


0.01* 


Peptoniphilus 


0.3 


(0-4.1) 


0 


(0-0.9) 


0.69 


Prevotella 


0.6 


(0-4.8) 


0 


(0-2.5) 


0.55 


Rhodanobacter 


0 


0 


0 


(0-0.9) 


0.16 


Sneathia 


0 


0 


0 


0 


0.07 


Staphylococcus 


0.1 


(0-0.3) 


0 


(0-0.3) 


0.34 


Streptococcus 


0.1 


(0-1.4) 


0 


(0-0.3) 


0.59 


Veillonella 


0 


0 


0 


(0-0.01) 


0.28 



FIG 3 Comparison of sequence-based urinary microbiome by cohort. The frequency (A) and median sequence abundance (B) of the overall most abundant taxa 
detected by sequencing were calculated. The families Enterobacteriaceae, Lachnospiraceae, and Ruminococcaceae could not be classified to the genus level. In panel 
A, a combination of Pearson chi-square and Fisher's exact tests was used to compare the frequency of genera detected by sequencing between the cohorts. 
represents a P value of<0.05.InpanelB,a Wilcoxon rank sum test was used to compare the median sequence abundances between the cohorts. IQR, interquartile 
range. "*" represents P values of <0.05. 
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TABLE 3 Richness and diversity measures of urinary microbiome 



UUI (n = 23)° 



Non-UUI control [n = 25)" 



Metric 



Total [mean (SD)] 



Mean (SD) 



Min-Max 



Mean (SD) 



Min-Max 



P value* 



No. of observed OTUS (richness) 97.02 (40.6) 86.35 (36.89) 25-155 106.84 (42.17) 39-201 0.08 

Chaol estimator (richness) 7,210 (5,331) 5,977 (4,647) 2,688-20,021 8,344 (5,750) 1,697-22,997 0.12 

Shannon index (diversity) 1.72 (0.96) 1.70 (0.91) 0.17-3.68 1.74 (1.02) 0.30-3.83 0.88 

Inverse Simpson index (diversity) 3.93 (4.28) 3.69 (3.39) 1.03-16.08 4.15 (5.07) 1.07-20.40 0.71 

17 Analysis performed on the subset of samples with detectable bacterial DNA. Min, minimum; Max, maximum. 
b Student's t test. 



these genera were isolated solely from the UUI cohort {Actino- 
baculum, Aerococcus, Arthrobacter, and Oligella); the concentra- 
tions ranged from 10 CFU/ml to 100,000 CFU/ml. Using matrix- 
assisted laser desorption ionization-time of flight mass 
spectrometry (MALDI-TOF MS), the majority of the cultured iso- 
lates were classified to the species level, one taxonomic level fur- 
ther than was achievable with V4 sequencing. This finer level of 
resolution revealed additional differences between the cohorts. 
For example, the genus Lactobacillus was frequently detected in 
both cohorts; however, when examination at the species level was 
done, distinct differences were detected between cohorts. Lacto- 
bacillus gasseri was more frequently isolated from the UUI cohort 
(27% UUI, 9% non-UUI; P = 0.02), whereas Lactobacillus crispa- 
tus was more frequently isolated from the control cohort (4% 
UUI, 18% non-UUI; P = 0.037) (Fig. 6). Also, Actinobaculum 
schaalii, Actinomyces neuii, Aerococcus urinae, Arthrobacter cum- 
minsii, Cory neb acterium coyleae, Gardnerella vaginalis, Oligella 
urethralis, and Streptococcus anginosus were more frequently iso- 
lated from UUI women (see Table S2 in the supplemental mate- 
rial). 



Comparison of culture versus sequencing for defining the 
urinary microbiome. More than half (58%; 30/52) of specimens 
examined by both EQUC and 16S rRNA gene sequencing tested 
positive for bacteria by both techniques. Bacteria were cultured 
from 14 (27%) sequence-negative urine samples, while bacterial 
DNA was sequenced from 3 (6%) EQUC-negative urine samples. 
Altogether, bacteria were detected in 90% (47/52) of urine sam- 
ples by EQUC or sequencing or both. Only 10% (5/52) of the urine 
samples tested negative for bacteria by both approaches (Fig. 7A) . 

Of the 30 urine samples that were positive by both sequencing 
and EQUC, there was considerable overlap in terms of the bacte- 
rial taxa detected (Fig. 7B). A total of 18 different genera were 
cultured, the majority of which were detected by sequencing, thus 
providing additional evidence of their presence in the urinary mi- 
crobiome and indicating that many of the sequenced genera rep- 
resent live bacteria. In some urine samples, EQUC detected genera 
that were not detected by sequencing; however, these nine genera 
were detected by sequencing in other urine samples, demonstrat- 
ing that these genera could be amplified by the universal primers. 
Trueperella was the only cultured genus that was not detected at 



UUI women 




I — I — I — I — I — I — I 1 — I — I — I — I — I 1 — I — I — I — I — I 1 — I — I — I — I — I 1 — I — I — I — I — I 1 — I — I — I — I — I — I - 
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Number of urine specimens 

FIG 4 Rarefaction curves of the cultured bacterial species by cohort. The plot depicts the number of species cultured via EQUC by the number of urine samples 
assayed. 
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FIG 5 Genus-level comparison of cultured urinary microbiota by cohort. The Pearson chi-square and Fisher's exact tests were used to compare the frequencies 
of the genera isolated from urine via EQUC. *, P <0.05; **, P < 0.001. 



the sequencing level in any of the urine samples. In contrast, some 
genera, such as Atopobium, were detected via sequencing but not 
by culture, suggesting that even the expanded culture technique 
was limited. 

DISCUSSION 

This analysis directly compared the urinary microbiome of 
women with UUI to that of women without UUI symptoms. We 
used two independent yet complementary techniques, 16S rRNA 
gene sequencing and extended culture, to characterize the female 
urinary microbiome. Both techniques identified statistically sig- 



nificant differences in the frequency and abundance of bacteria. 
These differences suggest a potential role for the urinary micro- 
biome in female urinary health and warrant further study. 

Each technique detected similar but not identical microbiome 
profiles. High-throughput sequencing provided the broader view 
of the bacteria present in the bladder, regardless of ability to cul- 
tivate the organisms. On the other hand, the cultured isolates 
could be identified to the species level, thus providing a finer level 
of resolution for the female urinary microbiome. Furthermore, 
these cultured isolates can be utilized in future studies to further 
investigate their potential for symbiosis or pathogenesis. Com- 
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bined, these complementary techniques yielded overlapping re- 
sults and provided the most complete description of the female 
urinary microbiome to date. 

In both cohorts, one or two genera dominated the majority of 
the sequence profiles. The two most frequently detected genera, by 
both sequencing and culture, were Lactobacillus and Gardnerella. 
The former are lactic acid-producing, facultative anaerobic bacte- 
ria known to play protective roles in the vaginal tract by decreasing 
pH and producing various bacteriostatic/cidal compounds (11, 
12). The latter are also facultative anaerobic bacteria frequently 
isolated in the vaginal tract, often in association with bacterial 
vaginosis. In this study, Lactobacillus was detected at similar fre- 
quencies in both cohorts but displayed lower median sequence 
abundance for women with UUI than for those without UUI. 
Furthermore, at the species level, the Lactobacillus species differed 



between cohorts. Whereas L. gasseri was more frequently cultured 
in samples from the UUI cohort, L. crispatus was more frequently 
cultured in samples from controls. Why the distribution of these 
two Lactobacillus species differs between women with and without 
UUI remains unknown, but it is possible that these two members 
of the normal vagina flora perform distinct functions in the blad- 
der. A similar scenario seems to be developing for Gardnerella, 
which was detected more frequently and in increased sequence 
abundance for the UUI cohort than for the non-UUI cohort. Since 
Gardnerella was detected in the urine of women without UUI, it is 
unlikely that the simple presence of Gardnerella indicates a dysbi- 
otic environment. G. vaginalis can be cultured from the vaginal 
tracts of women with and without bacterial vaginosis (13, 14). 
G. vaginalis strains isolated from these two groups of women vary 
in their ability to adhere and induce cytotoxicity, suggesting the 
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possibility of pathogenic and symbiotic strains of G. vaginalis ( 14) . 
Whether the isolates identified by MALDI-TOF MS as Gardnerella 
species fall into the latter category remains to be determined. 

Several bacterial genera were more frequently sequenced and 
cultured from the urine of women with UUI, including Actino- 
baculum, Actinomyces, Aerococcus, Arthrobacter, and Oligella. In- 
terestingly, many of these genera contain emerging uropathogens, 
including Actinobaculum schaalii, Aerococcus urinae, Oligella ureo- 
lytica, and Arthrobacter cumminsii (15-19). Whether these bacte- 
ria contribute to UUI is unknown at this time. 

In addition to Lactobacillus and Gardnerella, many of the gen- 
era detected in the urinary microbiome are often found in the 
vaginal tract, including Bifidobacterium, Enterococcus, Actinomy- 
ces, Prevotella, and Atopobium (13, 20). Since the urine samples in 
this study were collected via transurethral catheter, as opposed to 
a voided approach, these genera are likely true inhabitants of the 
urinary tract and not vaginal contamination. This contention is 
supported by our previous demonstration that the microbiome 
sequence profile of catheterized urine is distinct from that of 
voided urine and instead closely resembles that of urine collected 
via suprapubic aspiration, which bypasses the vagina and urethra 
(1). Genera found in both the urinary and vaginal tracts could 
suggest a shared urogenital microbiome. Alternatively, there 
could be differences at the species or strain levels, such that the 
vaginal and urinary community members differ. It is easy to imag- 
ine that the specific conditions of these quite different environ- 
mental niches could select for different traits. 

Two genera, Cory neb acterium and Streptococcus, were cultured 
more frequently from UUI women yet were detected at similar 
frequencies by sequencing. A third genus, Staphylococcus, was 
more frequently cultured from UUI women but more frequently 
detected by sequencing in the control cohort. The simplest expla- 
nation for these discrepancies would be that certain species of 
these genera do not grow under the current EQUC conditions. 
Efforts to further extend this new protocol are planned. 

Although the urinary microbiome shares a number of genera 
in common with other microbial sites within the human body, 
two key differences set the urinary microbiome apart. Compared 
to the colon, which contains 10 11 to 10 12 CFU/g (21), on average, 
the urinary tract microbiome consists of < 10 4 CFU/ml total bac- 
teria. In our study, the median amount of bacteria per urine was 
85 CFU/ml. The urinary microbiome, with a median inverse 
Simpson index of 2.3 and Shannon index of 1.5, is also less diverse 
than other microbial sites in the human body, such as the skin, 
mouth, and gastrointestinal tract (22, 23). Taken together, these 
findings suggest that the bladder is a relatively unique microbial 
site within the human body and may be more akin to other low- 
abundance sites, such as the eye (24). 

Changes in microbial diversity within a niche have been linked 
to disease. Examples include decreased fecal microbiome diversity 
associated with Clostridium difficile infection (25) and increased 
vaginal microbiome diversity associated with bacterial vaginosis 
(26). In this study, we did not detect large differences in sequence- 
based diversity between the two cohorts. It is possible that we were 
unable to detect a difference due to the lack of power with our 
small sample size. Alternatively, the amount of diversity within the 
urinary microbiome of UUI women may not differ from that for 
non-UUI women. Instead, the key distinction might be the bacte- 
rial types that are present. Although we did not detect differences 
in richness or diversity by sequence, we detected differences in the 



species richness between the cohorts by culture. A greater number 
of bacterial species were cultured from the UUI cohort than from 
the non-UUI control cohort. This finding further demonstrates 
the power of utilizing multiple approaches to define a microbial 
community. 

The urine of UUI women was more likely to contain Actino- 
myces, Aerococcus, and Gardnerella and less likely to contain Lac- 
tobacillus than urine collected from women without UUI. The 
mechanism behind this observation is unknown. One possibility 
is that the UUI bladder selects for some bacteria over others, and 
as such, the presence of these organisms in the bladder could serve 
as a marker for dysbiosis. Another possibility is that these bacteria 
contribute to UUI symptoms, a supposition supported by the ob- 
servation that each of the genera associated with the UUI cohort 
contains at least one reported pathogenic species. Taken together, 
this study, along with others, will eventually allow us to define a 
core or common urinary microbiome that can be used to detect 
alterations to that community. 

A limitation of this analysis is that our UUI and non-UUI co- 
horts differed in several characteristics that may have clinical rel- 
evance for the female urinary microbiome. Thus, we are unable to 
say whether age, BMI, or hormonal status affects our findings, and 
it is possible that our findings are related to these differences 
rather than to urinary symptomatology. Studies that describe the 
female urinary microbiome in large, well-characterized popula- 
tions have not been published. Future studies, with larger sample 
sizes, will be required to evaluate these potentially important dif- 
ferences. Studies describing the longitudinal stability of the female 
urinary microbiome (including populations undergoing treat- 
ment for urinary disorders) are also lacking. Such studies should 
be prioritized given the emerging evidence that the female urinary 
microbiome may yield important clinical information. 

MATERIALS AND METHODS 

Study design and population. Following Loyola institutional review 
board (IRB) approval, participants gave verbal and written consent for 
chart abstraction and urine collection with analysis for research purposes. 
Participants were recruited from the clinical practice of the Female Pelvic 
Medicine and Reconstructive Surgery Center of Loyola University Medi- 
cal Center between August 2012 and February 2014. These included 
women undergoing UUI treatment (UUI cohort) and a comparison 
group of women not bothered by urinary symptoms (non-UUI control 
cohort). All women were screened for potential study participation using 
the validated symptom questionnaire, the Pelvic Floor Distress Inventory 
(PFDI) (27, 28). Exclusion criteria for both cohorts included current UTI 
(based on urine dipstick) or history of recurrent UTI, antibiotic exposure 
in the past 4 weeks for any reason, immunologic deficiency, neurological 
disease known to affect the lower urinary tract, pelvic malignancy or ra- 
diation, untreated symptomatic pelvic organ prolapse (POP) greater than 
POP-Q stage II (vaginal protrusion more than 1 cm outside of the vaginal 
hymen), or pregnancy. Clinical and demographic information were ab- 
stracted from the electronic medical record. Enrolled participants com- 
pleted the long form of the PFDI and the Overactive Bladder Question- 
naire (OAB-q) (29). Premenopausal woman and postmenopausal women 
taking any form of estrogen replacement were considered estrogen posi- 
tive. 

Sample collection. Urine was collected aseptically via transurethral 
catheter. A portion of the urine sample was placed in a BD Vacutainer Plus 
C&S preservative tube for culturing. A separate portion for sequencing 
was placed at 4° C for <4 h following collection; 10% AssayAssure 
(Thermo Scientific; Waltham, MA) was added before storage at —80° C. 

The first urine samples collected (14 UUI samples and 8 non-UUI 
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control samples ) were screened exclusively by 1 6S rRNA gene sequencing; 
this was before we implemented EQUC. Subsequent urine samples (22 
UUI and 30 non-UUI control samples) were screened by both 16S rRNA 
gene sequencing and EQUC. The final urine samples (23 UUI and 15 
non-UUI control samples) were screened only by EQUC without high- 
throughput sequencing. Thus, we used EQUC for 45 UUI and 45 non- 
UUI samples and 16S rRNA gene sequencing for 36 UUI and 38 non-UUI 
samples. 

Urine culture. For standard urine culture, we struck 0.001 ml of urine 
onto 5% sheep blood (BAP) and MacConkey agars (BD BBL prepared 
plated media), which were incubated aerobically at 35°C for 24 h. Each 
separate morphological colony type was counted and identified in any 
amount. The detection level was 1,000 CFU/ml, represented by 1 colony 
of growth on either plate. If no growth was observed, the culture was 
reported as "no growth," indicating no growth of bacteria at the lowest 
dilution, i.e., 1:1,000. 

For EQUC, we struck 0. 1 ml of urine onto BAP, chocolate and colistin, 
naladixic acid (CNA) agars (BD BBL prepared plated media), which were 
incubated in 5% C0 2 at 35°C for 48 h (see Fig. SI in the supplemental 
material). A second set of BAPs was inoculated with 0.1 ml of urine and 
incubated in room atmosphere at 35°C and 30°C for 48 h. We also inoc- 
ulated 0.1 ml of urine onto each of two CDC anaerobe 5% sheep blood 
agar plates (BD BBL prepared plated media) and incubated either in a 
Campy gas mixture (5% 0 2 , 10% C0 2 , and 85% N) or under anaerobic 
conditions at 35°C for 48 h. The detection level was 10 CFU/ml, repre- 
sented by 1 colony of growth on any of the plates. Finally, to detect any 
bacterial species that may be present at quantities lower than 10 CFU/ml, 
1.0 ml of urine was placed in thioglycolate medium (BD BBL prepared 
tubed media) and incubated aerobically at 35°C for 5 days. If growth was 
visually detected in the thioglycolate medium, the medium was mixed, 
and a few drops were plated on BAP and CDC Anaerobe 5% sheep blood 
agars for isolation and incubated aerobically and anaerobically at 35°C for 
48 h. Each morphologically distinct colony type was isolated on a different 
plate of the same medium to prepare a pure culture that was used for 
identification. Matrix-assisted laser desorption ionization-time of flight 
mass spectrophotometry (MALDI-TOF MS) with the MALDI Biotyper 
3.0 software program (Bruker Daltonics, Billerica, MA) was used to iden- 
tify the bacterial isolates, as described elsewhere (4). To determine the 
false-negative rate, the following equation was used: the number of false 
negatives (EQUC positive, standard negative) divided by the sum of the 
number of true positives (EQUC positive) and false negatives (EQUC 
positive, routine negative). 

DNA isolation from urine. We used a previously validated DNA ex- 
traction protocol developed for the Human Microbiome Project. The 
protocol includes the addition of the peptidoglycan-degrading enzymes 
mutanolysin and lysozyme, which ensure robust lysis of Gram-positive 
and Gram-negative species, to isolate genomic DNA from urine samples 
(30). Briefly, 1 ml of urine was centrifuged at 13,500 rpm for 10 min, and 
the resulting pellet was resuspended in 200 fx\ of filter-sterilized buffer 
consisting of 20 mM Tris-Cl (pH 8), 2 mM EDTA, 1.2% Triton X-100, 
and 20 /ng/ml lysozyme and supplemented with 30 pi of filter-sterilized 
mutanolysin (5,000 U/ml; Sigma-Aldrich, St. Louis, MO). The mixture 
was incubated for 1 h at 37°C, and the lysates were processed through the 
DNeasy blood and tissue kit (Qiagen, Valencia, CA) according to the 
manufacturer's protocol. The DNA was eluted into 50 fx\ of buffer AE, 
pH 8.0, and stored at -20° C. 

16S rRNA gene library generation and MiSeq sequencing. Sequenc- 
ing was performed using a MiSeq desktop sequencer (Illumina, San Di- 
ego, CA). First, a 16S rRNA gene amplicon library was generated via two 
consecutive PCR amplifications. In the first reaction, the variable 4 region 
(V4) ofthe 1 6S rRNA gene was amplified using the universal primers 51 5F 
and 806R, which were modified to encode the Illumina MiSeq sequencing 
primer sequence at the 5 ' end (see Table S 1 in the supplemental material) . 
Reaction mixtures were incubated at 94°C for 2 min to denature the DNA 
template and amplified for 30 cycles at 94°C for 30 s, 55°C for 30 s, and 



72°C for 90 s. To ensure complete amplification, samples were incubated 
at 72°C for an additional 10 min. Ten-microliter aliquots of each reaction 
mixture were run on a 1% agarose gel. Samples containing a band of 
approximately 360 bp were considered PCR positive and subjected to 
further library preparation. Samples with no visible amplified product 
were considered PCR negative and not processed further. The PCR- 
positive reaction mixtures were diluted 1:50 and amplified for an addi- 
tional 10 cycles, utilizing primers encoding the required adapter se- 
quences for Illumina MiSeq sequencing and an 8-nucleotide (nt) sample 
index (see Table SI), using the PCR conditions described above. Unincor- 
porated nucleotides and remaining primers were removed via use of the 
QIAquick PCR purification kit (Qiagen, Valencia, CA), and the DNA 
concentration of the eluted product was determined by Nanodrop spec- 
troscopy (Thermo Scientific; Waltham, MA). One hundred nanograms of 
each sample amplicon was pooled and run through a 1 % agarose gel. The 
final product, which includes the V4 region and the adapter sequences 
(390 to 450 bp in length depending upon the length of the V4 region), was 
gel extracted via the QIAquick gel extraction kit and further purified via 
Agencourt AMPure XP-PCR magnetic beads (Beckman Coulter, Pasa- 
dena, CA) . The final concentration of the pooled DNA was determined via 
Nanodrop spectroscopy and diluted in EBT (Elution Buffer with Tris; 
Illumina, San Diego, CA) to 2 nM. An equal volume of 0.2 N NaOH was 
added, incubated at room temperature for 5 min, and quenched with 
hybridization buffer (Illumina, San Diego, CA) to a final concentration of 
8 pM. The Human Microbiome Project mock community HM-782D 
(BEI Resources, ATCC, Manassas, VA), a standard used to optimize our 
approach, was mixed 1 : 1 with randomly generated PhiX libraries, which 
are added to help focus the cameras on the sequencing clusters. This 
mixture was added to the sample library at equal volumes and placed in 
the 2 X 250 bp sequencing reagent cartridge according to the manufac- 
turer's instructions (Illumina, San Diego, CA). 

Care was taken to avoid bacterial DNA contamination by utilizing 
DNA-free reagents when applicable, filter sterilizing all solutions through 
a 0.2 ^M filter, and working in a PCR-clean hood. To control for the 
introduction of contaminating DNA, negative controls for extraction (no 
urine) and PCR (no template) were included in each experiment. The 
extraction negative control for each experiment was sequenced to identify 
spurious genera likely introduced from contaminated reagents and mate- 
rials. 

Sequence analysis. The Illumina MiSeq postsequencing software pre- 
processed sequences by removing primers and sequence adaptors. Using 
the open-source software program mothur (v 1.31.2), the paired-end 
reads were assembled and contigs of incorrect length (<285 bp or 
>300 bp) and/or contigs containing ambiguous bases were removed (31, 
32). These modified sequences were aligned to the SILVA reference data- 
base, and any potential chimeric sequences were detected and removed 
using the program UCHIME (33). The remaining sequences were taxo- 
nomically classified, using a naive Bayesian classifier (34) and the mothur- 
formatted RDP training set v9, and clustered into operational taxonomic 
units (OTUs), an operational definition of a species-level cluster based on 
sequence similarity, using a 97% cutoff. The software program METAGE- 
Nassist was used to link OTU nomenclature to taxonomic assignments 
(35). 

All samples were processed in duplicate, and the percent reads of the 
replicates were averaged for downstream analysis. Stacked bar plots based 
on sequence abundance were produced for each sample. Euclidean dis- 
tance was calculated between samples, and the complete method was used 
for hierarchical clustering via R software, version 2.15.1 (36). Richness 
and diversity metrics, including the number of observed OTUs, Chaol 
estimator, Shannon index, and inverse Simpson's index, were calculated 
using mothur and were based on subsampling to the number of sequences 
in the sample with the least coverage. Urine is a low-biomass environ- 
ment, and thus the sequencing results are more likely to be influenced by 
extraneous DNA arising from reagents and resources; therefore, only 
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reads representing >0.01% of the sample total were included in the 
stacked bar plots, dendrogram, and frequency and abundance analyses. 

Statistical analysis. Statistical analyses, comparing participant demo- 
graphics and symptoms, were performed using the SPSS software pro- 
gram, version 19. For continuous variables, Student's t tests were applied. 
For categorical variables, Pearson chi-square and Fisher's exact tests were 
performed. Results were considered significant when the P value was less 
than 0.05. 

Statistical analyses of the microbiome data were performed using the 
SAS software program, version 9.3. The Wilcoxon rank sum tests were 
used to compare the median abundances for the 15 most abundant se- 
quenced taxa and all cultured genera between groups. The frequencies of 
detected genera were compared between groups, using either the Pearson 
chi-square or Fisher's exact test, depending on assumption validity. No 
adjustments for multiple comparisons were made, since these analyses 
were considered descriptive. 

SUPPLEMENTAL MATERIAL 
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